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We numerically investigate the decay, via parametric res- 
onance, of the inflaton with an nn?(fp' potential into a scalar 
matter field with a symmetry breaking potential. We con- 
sider the case where symmetry breaking takes place during 
inflation. We show that when expansion is not taken into ac- 
count symmetry restoration and non-thermal defect produc- 
tion during reheating is possible. However in an expanding 
universe the fields do not spend sufficient time in the instabil- 
ity bands; thus symmetry restoration and subsequent domain 
wall production do not occur. 
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Introduction 

The realisation that the decay of the inflaton might 
occur explosively, during a stage dubbed "pre-heating" 
10-^ and lead to a universe far from thermal equilib- 
rium, has a number of important ramifications. The most 
obvious is that the temperature of the universe, after 
thcrmalisation, may be much higher than under the old 
model of reheating. Another is the recent construction 
of a model for baryogenesis Q which takes advantage of 
the out-of-equilibrium nature of the universe shortly af- 
ter pre-heating. We, on the other hand, will be interested 
in the possibility of non-thermal phase transitions after 
pre-heating, and the possible production of topological 
defects. 

This suggestion was first made by Kofman et al. 
and Tkachev ||]. The essential idea is that pre-heating, 
via the mechanism of parametric resonance, typically 
gives rise to very large fluctuations in the fields and it is 
these that may restore a symmetry which has been bro- 
ken by the time inflation ends. The important point is 
that symmetry restoration may occur even if the symme- 
try breaking scale is higher than the final reheat temper- 
ature. Then as the universe subsequently expands and 
cools, phase transitions will result in the creation of topo- 
logical defects. This would be very intriguing. For a start 
we might find the re-emergence of a problem inflation 
was designed to solve, namely the monopole problem. 
The creation of domain walls would also be potentially 
problematic. However the production of strings may be 
a desirable feature of such scenarios, in that strings can 
be a seed for structure formation. 



Khlebnikov et al. |7| have recently demonstrated that 
non-thermal defect production indeed occurs during in- 
flationary reheating. They considered a model in which 
the inflaton field has a double well potential and is 
coupled to massless scalar matter fields Xi,i = 1 . . .N. 
Similar results (obtained from 2-D simulations) have also 
been shown by Kasuya and Kawasaki |^ . 

We study a different model. In our work the inflaton 
has a unique minimum and decays via parametric reso- 
nance into a scalar matter field x with the discrete sym- 
metry X —X- That is, we are interested in the question 
as to whether domain walls are produced. (Note that we 
are not considering a two-field model of inflation.) We 
assume that symmetry breaking occurs during inflation, 
for if it occurs afterwards the symmetry will then be re- 
stored a number of times by the oscillation of the infla- 
ton alone, without any need for parametric resonance, 
and defect production will almost certainly ensue. This 
was pointed out by Kofman et al. in their original paper 
Q, where this model was first examined in the context 
of parametric resonance and defect production, and by 
Kofman and Linde in [^. 

We investigate both the expanding and non-expanding 
cases and we show that symmetry restoration and defect 
production only occur in the non-expanding case. How- 
ever the defects are not stable. Our work is based on 
3-D lattice simulations, though we also understand the 
results from an analytical point of view. 

The paper is set out as follows. We begin with a de- 
tailed explanation of the model including necessary con- 
straints on the model. The second section is devoted to 
the analytical understanding of the model and predic- 
tions for the full simulation. After this we highlight the 
workings of the numerical code and present the results 
from our simulations. The final section is a brief sum- 
mary of our study and conclusion. 

1. Model 

We consider the simplest model of infiation and reheat- 
ing: the inflaton </) has potential im^t^^ and it decays into 
a second scalar field x- We let x have the simplest po- 
tential which allows symmetry breaking: \\{x^ ~ Xo)^' 
and suppose it interacts with through the interaction 
term: \g'^(l?"X^ ■ Thus the full potential for the two-field 
theory is: 



X^)^ + ^9'<^\' (1) 
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The most important constraint on our model is the 
requirement that the symmetry breaking occurs during 
inflation. This is equivalent to demanding the interaction 
term in (|l|) not lift the degeneracy of the potential in 
the X direction. As (f) is homogeneous to a very good 
approximation for initial times, we may write (j) ~ (j){t) 
so that the effective potential for x is: 

V^"^lx{x'~xir + l9'^\t)x' (2) 

This has degenerate minima if Axo > g'^'Pit)'^- If this is 
true at t = 0, then it is always true; for in the expanding 
case (f) is redshifted to zero and in the non-expanding case 
(j> is not observed to increase. If we define (po = (j>{0), then 
the condition for degenerate minima becomes: 

Axg > g'^l (3) 

Another condition follows from the requirement that 
inflation proceed as usual, despite the presence of x- For 
this we must have the vacuum energy of x to be much 
smaller than that of 0: 

\xxt> « Irn^^l (4) 
When is combined with (^), we obtain: 

5^X0 < (5) 

which ensures the frequency of oscillation of (f> after in- 
flation is dominated by m. 

We next demand that the maximum reheat tempera- 
ture be insufficient to lead to a thermally-corrected po- 
tential for X with the degeneracy lifted. Since Trjj < 
(m^o)"^! this condition becomes: 

Xo>( — )^9!'o (6) 

Actually the final reheat temperature will be many orders 
of magnitude smaller than we have indicated, so (^) is 
much stricter than it needs to be. However it is easy 
to satisfy and we will retain it in order to stress we are 
interested in non-thermal defect production. 

Finally there are the usual constraints on m and 
arising from the COBE data: m < lO^^Mp and g"^ < 
10~^; the latter ensures radiative corrections do not lead 
to a too large self-coupling for (though in a super- 
symmetric model of inflation this condition may not be 
necessary ]lO|). 

2. Analysis 

Our working hypothesis is that parametric resonance 
will generate large fluctuations in x which will restore the 
symmetry of the x-field in some regions of space. After 
this, in such regions, x may evolve into the two degen- 
erate minima of its potential, possibly forming a stable 



domain wall. In other words, we expect the mechanism 
of parametric resonance will be able to channel sufficient 
energy from the infiaton into the x-field in order to drive 
X over the potential barrier. In this section we outline the 
conditions for effective parametric resonance, considering 
both the expanding and non-expanding cases. 

We begin by assuming the usual flat FRW universe, in 
co-moving co-ordinates: 

ds"^ = dt'^ - a^{t){dx^ + di/ + dz^) 

To a good approximation the coherent oscillations of 
the inflaton rapidly give rise to a matter dominated uni- 
verse so we may write: a{t) = (1 + •mt)'s ^ where i = is 
the time at the end of inflation and we have used our free- 
dom to scale the spatial co-ordinates to impose a(0) = 1. 
If we put: a{t) = (1 + mt)", we can simultaneously con- 
sider the non-expanding case [n = 0) and the expanding 
case {n = 

The equations of motion for the flelds and x are the 
standard ones for scalar fields in a curved space-time with 
minimal coupling to gravity: 

4> + 3H^-a-'^V^(j)+{m^ +g^X^)(l)^0 (J) 
X + 3HX- a-^V\ + (A(x' - xl) + 9^4>^)x = (8) 

where H = -. 

a 

These are the equations which are solved in our numer- 
ical simulation. Here we point out that in our simulation 
the scale factor a{t) is put in by hand, rather than deter- 
mined in a self-consistent manner. 

The initial conditions for <j) are those of the end of in- 
flation i.e. when the slow roll approximation is no longer 
valid. We find: ()i(0) = (pa ^ vfe^'^^ "^C^) = -mcpo. 
In addition we chose m = lO~^0o- For the non-expanding 
case we will assume 0(0) = 0. As for X: inflation dic- 
tates it will be initially comprised of quantum fluctua- 
tions about one of the instantaneous minima of V-^^-^, 

which are given by x±(0 = ^Xo Vi^~ x% ) ■ Without 
loss of generality we may choose the positive minimum. 

To make headway in our analytical understanding of 
the model we must make some simplifying approxima- 
tions. As already noted, we may take to be homoge- 
neous to start with. Recalling (||), this allows us to solve 
(^): = 00 a,{t)^2 cosmt, and to re- write (^) as: 

X + 3Hx-a-^V\+-^V^" = (9) 

In order to understand the evolution of the x-field, at 
least initially, we will expand V^-^-^ to quadratic order 
about x+{t)- 



^The results of this section are valid only for these two values 
of n. 
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xXit)) 

(10) 



How now do we ask if symmetry restoration occurs? 
We note that U^{^x+{t)) = j^Xo> which is the barrier 
height of V^^^ . Thus the question as to whether x sur- 
mounts the potential barrier of V^-^^ may be replaced 
with the question: does Xi evolving under the influence 
of C/^, ever become less than i^+C*)? 

With the substitution of V^^^ with U-^ in (^), we ob- 
tain a linear equation for x- Upon rescaling: {t — > 
X ^ ^, ^ X - <^0X, 9^ - {f,?9\ A ^ 

(^)^A}, taking the Fourier transform and introducing 

Xk = a5 Xkj we arrive at a Mathieu equation with (in 
general) time-dependent co-efficients, for the non-zero 
momentum modes X\^: 



+ [Ak {t) - 2q{t) cos 2i) Xk = 



(11) 



where = ^ + 2\xl - ^ andq^ 

The existence of exponentially growing, or unstable 
solutions to the Mathieu equation is well knowrj^ and, 
in physics literature, has been termed parametric res- 
onance. However recent work [p^|-p^, inspired by the 
pre-heating scenario, has led to new insight into the na- 
ture of these solutions. For time-independent coefficients, 
we may identify three different regions of instability, (i) 
Narrow band resonance: Ak P > 2q, I = 1,2... and 
q <C 1. The solution Xk oscillates with period 2tt, and is 
modulated by an exponential function, e'^**, where is 
called the characteristic exponent or growth index. Typ- 
ically ~ |. (ii) Broad band resonance: Ak > 2q and 
g > 1. The solution oscillates many times faster than 
in the narrow band case and grows in exponential jumps, 
occurring at time intervals of tt. The overall growth index 
may also be much larger: one finds |/ife| < 0.3. (Hi) Neg- 
ative instability: Ak < 2q. The solutions in this region 
are almost all unstable and > 1 is possible. Recently 
Greene et al. have constructed a model which makes 
use of this region. 

The situation of the Mathieu equation with time- 
dependent coefficients is much more complicated. The 
difficulty is that, during one oscillation of the driv- 
ing force, Xk does not remain in one instability band 
but in fact passes in and out of many such bands. 
The result is that parametric resonance becomes essen- 
tially haphazard — may even decrease in amplitude 
at times — and hence the term stochastic resonance has 
been applied to this scenario. However exponential type 
growth is still possible, and the work of Kofman et al. 
[Of has done much to clarify this issue. 



Returning to equation ([1 we see that in our case Ak 
is always greater than q because of (||) , hence we will only 
see broad band or narrow band resonance. A necessary 
condition for parametric resonanc efl @ in these regions 
is: 



^2 _ 



Ak 



2V9 



< 1 



(12) 



It is useful to introduce the dimensionless parameters: 
a = g'^, (3 = Xxl and 7 = Xo^- Then conditions (||), (Q) 
and (^) may be succinctly written: a < ^ 7 < 10^. 

For n = |, k2 ^ 



t which must always become 



larger than 1. Later we will see that the typical time for 
symmetry restoration to begin in the non-expanding case 
is i 30. Thus to satisfy ( p^ ) in the expanding case we 
will require ^/a 3> (3. But a < f3. So we must have: 



\/a < 1 



(13) 



Or in other words, gp^yg ^ 10^^^. We conclude that if 
disl ) is not satisfied, and we believe it is not in any re- 
alistic model, then the system is not long enough in the 
instability regions to restore the symmetry. The conclu- 
sion is further reinforced if we realise we also require^ 
fik > ^H, but that initially |i?(0) = 1 > 0.3 ~ f^k,max- 

The necessary condition for parametric resonance in 
the case of n = is much less stringent. (O) gives: 



(14) 



A test of the validity of our approximation in the non- 
expanding case is the time taken for symmetry to start 
being restored. As mentioned previously, what we need 
to calculate is the time at which x first becomes less than 
^X+{t)- The modes in the resonance band will dominate 
the behaviour of x so: 
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X 



(27r) 2 (1)0 J res. band 



d^kx^e 



-i k.x 
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Afc 



(27r)t (/-o 



'res 2 IXk, res I 



(27r)2 00 y/LOres 

where Afc is the width of the resonance band and Wres = 
Ak,res — 2gcos2f = fc^gj + 2f] — a — acos2t. The ap- 
pearance of ^ and y/ ujres is due to x starting out as 
quantum fluctuations ||l5[| . 



^For a good review see p 



■^This is a more generous condition than Eq. (56) in [ Il3[ , 
where 7r~^ replaces 1 on the RHS. We have found, for our 
parameter choices, that (^2|) includes the (two) main broad 
band instability regions. 

*This follows from a simple-minded calculation of the ratio 
of energy gained by Xk to energy lost due to expansion. 



3 



Typically in our parameter range, we find a region of 
broad band resonance for [3 slightly greater than a, with 
kres ~ Afc ~ 1 and ^ 0.25. Then the statement 



X < \x+{t) becomes ^ f]' 



Xo, which implies: 



t>4:\n{10^(3ij-2) - 30 



(15) 



To summarise: our analysis suggests that in the non- 
expanding case the conditions for parametric resonance 
depend only on a and whereas the time taken for 
parametric resonance to restore the symmetry of the x~ 
field is determined mainly by 7. In the expanding case 
we do not expect domain walls to be produced. We now 
turn to the full simulation of the two-field system. 

3. Simulation and Results 

The classical equations of motion (0) and (^) are a good 
approximation to the behaviour of an actual matter field, 
X, coupled to the inflaton, provided the mode amplitudes 
of X are large. Although initially these are taken to be 
quantal fluctuations, the modes we are interested in are 
in the resonance band and hence they quickly grow and 
become classical |^^. 

We work in terms of the rescaled, non-dimensional 
quantities of the previous section. The universe becomes 
the usual 3-dimensional lattice in co-moving co-ordinate 
space, and we apply periodic boundary conditions. The 
fields are evolved using a leap-frog scheme, which is an 
explicit algorithm second-order accurate in time. The 
spatial derivatives are evaluated to fourth-order accuracy. 

The results presented here were obtained on 64^ lat- 
tices, though we also checked representative cases on 128"^ 
lattices. We varied the lattice spacing, the only pro- 
viso being that we remained sensitive to the momentum 
modes in the expected resonance bands, and obtained es- 
sentially identical data. The total energy of the system 
was conserved to better than 0.1%. Two other tests con- 
vinced us of the accuracy of our simulation: we were able 
to repeat the results of [0, and, in an 0{2) x-ficld sim- 
ulation, we found a faithfully modelled, collapsing string 
loop. 

A brief note about the determination of the initial con- 
ditions for X is warranted. As mentioned before, we treat 
X initially as quantum fluctuations about one of the in- 
stantaneous minima of its potential. One may think of 
the quantum fluctuations in momentum space as an infi- 
nite collection of harmonic oscillators (with time depen- 
dent frequencies in this case) ||l^. Our approach essen- 
tially was to pick the occupation numbers of the modes on 
the lattice randomly from a Gaussian distribution with k- 
dependent width |l^ . We also included a random phase 
and then inverse Fourier transformed to obtain x in con- 
figuration space. We checked that our results were not 
sensitively dependent on the particular random numbers 
chosen. A final important point is that the lattice spac- 
ing gives a natural cutoff to the high momentum modes. 
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FIG. 1. The mean and variance of x in time (q = 100, 
(3 = 130 and 7 = 10''). 



and one must choose this so that Efiuct. ^ Einfiaton- 

We examined the region in parameter space of a, ~ 
100 and 7 ^ 10**; in other words g'^f^y ~ lQ~^^,Xphy ~ 
10-6 and xo^phy ^ 0.01(/)o ~ 3 x lO^^ GeV. The choice 
of these parameters was partially determined by the re- 
quirements that we be sensitive to the instability bands 
and the box be large enough to include any defects pro- 
duced |po| . The numerical parameters were: lattice spac- 
ing Aa; = 0.1 — 1 and time step Ai = 0.02. 

The first result is that we did not find any evidence 
of symmetry restoration in an expanding universe. In all 
cases we considered, x ~* Xo and (j) merely redshifted to 
0. This is as predicted by ( p^ and Kofman et al. ||^. 

The non-expanding scenario is more interesting how- 
ever. Figure 1 shows the behaviour of (x) and ^/((x^) ~ 
(x)^) for a = 100, /3 = 130 and 7 = 10^. This is typical of 
the cases in which symmetry is restored: the mean value 
of X ^ and the fluctuations grow to be of the order of 
Xo- Also typical is the fact that the back-reaction of x 
on (/) is very slight. (^) continues to oscillate sinusoidally 
and the fluctuations of are only of order 3 x lO^"^ <C 1 
at t = 200. 

Plotted in flgure 2 is the fraction, /_ , of the volume of 
the box with x < 0. The jumps in /_ occur at intervals of 
TT which is as expected; tt is the period of the driving term 
in the Mathieu equation i.e. c\? . Note too the onset of 
symmetry restoration is at i « 30, conflrming our earlier 
order of magnitude estimate. However, considering all 
our runs, we did not flnd agreement with the speciflc form 
of (15). (Usually t — 15 — 80.) On the other hand varying 
7, for fixed a and did not alter whether symmetry was 
or was not restored, as was our prediction. 

In the non-expanding case, symmetry restoration 
(when it occurs) gives rise to non-thermal production 
of defects. However these defects are not stable. This 
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FIG. 2. The fraction of the volume with x < (a = 100, 
^ = 130 and 7 = 10^). 
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FIG. 3. The mean of x in time {a = 100, /3 = 110 and 
7 = 10^). 

can be seen in our simulations: the regions of x < are 
essentially randomly scattered throughout the box (even 
when they account for half the volume of the box) and 
percolate. In addition they are uncorrelated in time. 

There is one problem with the particular case we have 
been considering above: it does not satisfy ([l^) . Thus the 
analysis of section 2 cannot be entirely correct. Figures 
3 and 4 detail what was an unexpected result but which 
led to a more complete understanding of our model. Here 
a = 100, 13 = 110 and 7 = 10"*. Clearly what is hap- 
pening is an initial series of global sign changes, prior 
to eventual local symmetry restoration. The change in 
sign is not important because the symmetry of x remains 
broken, but it does indicate we were naive to imagine x 
would fluctuate about ^'^^ early times. 

It is more appropriate to put x(x, t) — Xh{t) + Sx{'x., t), 
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FIG. 4. The fraction of the volume with x < (a = 100, 
13 = 110 and 7 = 10''). 

where Xh satisfies the homogeneous version of §). Then 
the linearised equation for X]^ = a^Sxk, which replaces 
(111), is: 

Ik + + H^xlit) - xg) + .g'</''W)^k = (16) 

The behaviour of Xh is quite delicate. Numerical so- 
lution of its ODE shows that for early times it is true 
Xh ~ X+7 except when /3 « a where we find x/i ~ ±X+- 
This is the sign change seen in figures 3 and 4. At later 
times Xh typically oscillates about x+, with a frequency 
several times greater than that of the inflaton. For n = |, 
Xh ~^ X+ ~^ Xo so in this case we believe our earlier 
analysis is essentially unchanged. However for n = 0, 
Xh tends to increase in amplitude. This ruins a straight- 
forward analysis in terms of a Mathieu equation for 
(which would come about if Xh ~ X+)- Numerical analy- 
sis of (^6|) for n ~ though does reveal exponential-type 
solutions. In fact parametric resonance appears to be 
much more likely: the instability bands are larger than 
those of the broad band regime of the Mathieu equation. 
We are currently working on an analytical understanding 
of these results Q . 

To conclude this section figure 5 shows the region of 
the a — f3 plane we investigated where symmetry restora- 
tion occurs, for the non-expanding case. We stress that 
this plot is an extrapolation from 45 simulations we per- 
formed, and as such to be taken only as a rough guide. 
However it indicates two features which we are certain 
about. Firstly the existence of band structure. This is 
reminiscent of the Mathieu equation instability bands, 
even down to their relative sizes. Secondly the fact that 
as P increases for fixed a, symmetry restoration is no 
longer possible. From this graph we conservatively de- 
duce a condition for symmetry restoration: 
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FIG. 5. The region of symmetry restoration (darker shad- 
ing) in the a — (3 plane. 



/3 < 2a 



(17) 



Also shown in figure 5 is our prediction from section 2. 
The discrepancy is the most striking evidence that our 
earUer analysis was inadequate. 

4. Conclusions 

We have numerically investigated the decay of the in- 
flaton into a scalar matter field x with a symmetry break- 
ing potential, when the symmetry is broken during infla- 
tion. In the non-expanding case we found that symme- 
try restoration is possible via parametric resonance. The 
analysis is complicated but interestingly the regions of 
instability are much larger than those of the Mathieu 
equation. This situation of more efficient parametric res- 
onance has also been seen in the work of Zanchin et al. 

We think it would be useful to more fully under- 
stand equations like (^|), perhaps utilising the methods 
of H. Finally we found that non-thermal defect pro- 
duction occurs, but that the defects are not stable. 

In the expanding case we found that domain wall pro- 
duction did not take place. This leads us to conclude 
that the inflaton may safely couple to x with the poten- 
tial as in (|l|). Actually the precise form for the potential 
may be unimportant. Recall that the symmetry of x is 
broken before the end of inflation. In fact, if there are 
to be no fluctuations due to the initial symmetry break- 
ing in our Hubble volume, the symmetry must be broken 
early enough to allow for the usual number of e-foldings 
by the time inflation is over. In our model this means 
/3 > 100 a [|0| which is to be compared with (|l^. This 
amounts (literally) to a very formidable barrier to sym- 
metry restoration. This suggests it is unlikely that there 
will be defect production in any realistic model in which 
X has the symmetry breaking potential and the symme- 
try is broken during inflation. 
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